Probing an untouchable environment as a resource for quantum computing 
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The full control of many-body quantum systems is expected to be a key towards the future nano- 
and quantum technologies. Among others, the realisation of quantum information processing jj has 
been studied intensively as a good test bed of quantum control as well as an ultimate engineering 
task that makes full use of quantum mechanical effects @, S|- Yet, manipulating quantum states is 
extremely hard, since information encoded in quantum states easily leaks out to the environment 
due to the inevitable interactions with it. In the theory of open quantum systems, an environment 
is usually treated as a large bath 0, [Hi , assuming that most of its dynamical details are averaged 
out. Thus we do not see it as a quantum object that we can control actively. In realistic situations, 
it is indeed near-impossible to precisely identify the quantum nature of an environment, let alone to 
control it at will. However, here we demonstrate how this formidable task can be achieved, provided 
the dimension of the environment can be regarded as finite within the timescale we manipulate 
the system. The information thereby obtained will be useful not only for deeper understanding of 
the system dynamics under decoherence but also for exploiting the environment as a resource for 
quantum engineering, such as quantum computation. 
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Specifically, our primary goal here is to identify the 
Haniiltonian Hse, which describes the effective interac- 
tion between the principal system S and the surrounding 
system E and (a part of) its internal dynamics. The di- 
mension ds of E is also unknown a priori and is a subject 
of our identification. We shall call E the environment 
symbolically throughout the paper. The knowledge of 
HsE is vital for controlling E a,s a, useful system, not to 
mention for checking its controllability. Readers may be 
reminded of the methods of quantum process tomography 
(QPT) [l|, 0-[ll| as means to determine all the parameters 
that characterise a general quantum evolution, namely a 
completely positive (CP) map. Nevertheless, QPT is a 
scheme to estimate the CP map for a quantum system 
for which we can prepare a specific state and perform 
measurements. Thus, the conventional QPT methods do 
not reveal the nature of an environment, which is beyond 
the reach of our measurement. 

Incidentally, there have been a series of studies on 
Hamiltonian identification of a many-body system un- 
der the condition of limited access |12h15| . However, all 
of them assume that a priori knowledge was available 
about the system configuration and the initialisability of 
the system state that depends on the type of interaction. 
In the present analysis, no particular assumptions are 
made about system structures or the type of interaction, 
except for the finiteness of ds- Therefore, the task is 
even more nontrivial than existing tomographic schemes. 

The basic idea is as follows. We first steer the joint 
system SE and an ancillary system A into a maximally 
entangled state with respect to the partition between 
SE and A, while S and E keep interacting through 
HsE- Then, the entire state on 5*, E, and A will be 
a product of two maximally entangled pairs of subsys- 
tems, which are S and Ai C A, and E and A2{— A\Ai). 
The time evolution of SE, expressed by an operator 
UsE{t) = e-xp(—iHsEt), can be observed as a mirror 



image on the side of Ai and A2, due to a property of 
maximally entangled pairs. The desired information on 
the Hamiltonian Hse can be extracted by performing 
QPT on S and A. 

We shall primarily focus on the ideal case where the 
effect of errors is negligible, in order to elucidate the 
essence of the idea. Particularly, the consideration of 
perfect entanglement between SE and A allows us to 
rigorously prove the propositions necessary for the main 
results. 

Let Hs and He the Hilbert spaces of the principal 
system (S) and its environment (E), whose dimensions 
are ds and d^, respectively. A key assumption we make 
here is that dE is finite, although it can be unknown. 
Naturally, dE niay be infinitely large in general, but we 
consider a situation where the system 5* effectively inter- 
acts with only a finite dimensional subspace E of the uni- 
verse E'. That is, the interaction between S and E is so 
dominant within the relevant timescale for describing the 
dynamics of the system that we can justify this assump- 
tion. In other words, the combined system in Hs ®He 
undergoes a unitary evolution. 

For longer timescales, the combined system SE cannot 
be immune to the effect of interactions with its surround- 
ing environment E' . A state psE on Hs ®He'^^ now 
subject to equilibration and will tend to some fixed state 
p^g^,. This fact can be used to reset the state psE before 
iterating the protocol. As will be clear later, the initial 
state p'^gl^ can be arbitrary, as long as it is a fixed state, 
albeit unknown. 

Another important assumption is that we are allowed 
to provide as many ancillary states as necessary, each of 
which is a maximally entangled pair between ai and a^'- 
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evolve for a time duration (< ilp ). 

FIG. 1: The protocol for steering the state on SEA towards 
the maximally entangled one between SE and A. The thick 
blue lines and the red dotted line represent entanglement and 
the interaction, respectively. The three systems, S, E, and 
A, are in some entangled state after foregoing rounds of the 
protocol, while j4 is a null space for the first round. In (a), a 
maximally entangled pair \Ta-^a2) is provided and the states 
of the subsystem ai and S are swapped to make the entan- 
glement network look like (b). Relabeling A, ai, and a2 as a 
new A as in (c), we perform state tomography of A and subse- 
quently a local filtering operation J-f*p on A. If J-f^p succeeds, 
iterate the procedure, going back to (a). 



They will form the ancillary system A as the protocol pro- 
ceeds. Also, the interaction between A and E is assumed 
to be negligible. Further, we shall take it for granted that 
PSE is initially pure, because, otherwise, we can always 
purify it by appending an additional Hilbert space. 

In order to make use of the 'mirroring effect ' of entan- 
glement to identify Hse, we first need to steer the state 
on SEA to establish maximal entanglement between SE 
and A. Let us describe how the state-steering protocol 
goes, and delineate why it works out for our purpose. 
Figure [1] depicts the state-steering protocol. We start 
with an initial (fixed, but unknown) state pg*^ and abun- 
dant copies of \Ta1a2) in Eq. (P). At t = 0, Ha is a 
null space, supporting no states. The SWAP operation 
between S and oi, which must be fast enough compared 
with the system dynamics, will be denoted as SWAPsai . 
The C-th round of the protocol proceeds as follows (C 
starts from zero): 

Step 1: Apply SWAPga^, where ai is a subsystem of the 
newly provided |Taia2)j ^^d then let A incorporate 
ai (the former S) and a2. 



Step 2: Perform state tomography of pA, and apply a local 
filtering operation J^^p on pA- If it fails, carry out 
the whole protocol from the beginning. 

Step 3: Let the SE system evolve for a time duration (< 
Ate) so that the functional of psA, ^Esa, which 
is defined below by Eq. ([2]), increases by ec > 0. 
See Sec. IV of the supplementary information (SI) 
as to how we should determine Ate and ec- 

Step 4: Terminate the protocol if AEsa is found to be non- 
increasing; otherwise, let the SE system evolve so 
that AEsA > ec, and go back to Step 1. 

Intuitively, Ate and ec are set so that we can com- 
plete the steering protocol within the desired time period, 
which can be made as short as possible at the expense of 
success probability. 

The local filtering operation on pA is written as 

^LF PA = ^lfPa^'lf: where Flf = \J Kiin ■ Pa^ with p'^'^ 
the inverse of pA on its support and Amin the smallest 
eigenvalue of pA- The probability of success of J^^p is 
Amin • rankpA- 

The quantity AEsa we measure in Step 3 is defined as 

AEsA := SipsA) - SipA) + Inds, (2) 

where S{p) = —TT{plnp) is the von Neumann entropy 
and Pa — TtsPsa- While what it represents may not be 
obvious at first sight, AEsa is the change in entangle- 
ment between SA and E due to SWAPsai in Step 1 of 
the following round. To see this, letting E{\ip)) be the 
amount of entanglement in state {ip) with the partition 
SE\A, we have 

AEsA = E{\^'sea'))~E{\'^sea)) 

= Eil'i'EAa,) <E> \Tsa2)) - EH'^'sEa) <E> |T,,a.)) 

= S{psa) - S{pa) +lrids, (3) 

where A' is the relabeled A, i.e., Aaia2. 
We can rewrite AEsa as 

AEsA = D{pse\\ps ® Pe) + D{ps\\I^lds), (4) 

where D{p\\a) = Tr(plnp— plnci) is the relative entropy 
between p and cr, and is a x ds identity matrix. 
Equation ([4]) shows AE > due to the nonnegativity 
of the relative entropy. Thus, we can see that, when 
AE is small, psE and ps are very close to the product 
state Ps ® Pe and the completely mixed state /ds, 
respectively. 

The above observation implies that when AEsa — 
there is a subsystem Ai of A that is maximally entangled 
with S. Further, we expect that the remaining part A2 
of A is maximally entangled with i? as a result of J-{^p. 
This naive guess is proven in detail in the SI. 

Note, however, that the argument in the SI involves 
some mathematical subtleties concerning multiple possi- 
bilities of the set (d^, \'^sea},E[se) that leads to identi- 
cal observable dynamics on S and A. We call this set of 
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FIG. 2: The desired state of S, E, and A. A subsystem Ai of 
A and its complement A2 are maximally entangled with 5* and 
E, respectively. Further executions of the entangling protocol 
do not change the entanglement structure: They only enlarge 
the size of A2 without increasing the entanglement with E. 

three ingredients a triple. More precise propositions we 
show in the SI are as follows: 

• There can be equivalent classes of triples, so that 
all triples within a single equivalent class give rise 
to the identical observable dynamics on S and A, 
regardless of any active controls on them. For our 
purpose of controlling E, it suffices to identify one 
in the class for the observed time evolution on SA. 

• The resulting state |^5_ba) of the state-steering 
protocol satisfies a condition, which is expressed 
as 

TiEil'i' sea) {"^seaD = \Tsa,){TsaA®Pa,, (5) 

where Ai and A2 are non-overlapping subspaces of 
A such that "Ha = Hai ® 'HA2, I'^SAi) is a maxi- 
mally entangled state, and is a projector onto 
'Ha2 ■ We refer to Eq. ^ as the maximal entangle- 
ment {ME) condition (Sec. IV of SI). 

• If the state I'^^sea) fulfills the ME condition, then 
there is a maximal entanglement between systems 
SE and A. That is, there exists an equivalent 
triple, in which A2 is fully entangled with the entire 
E. Thus 

I'^sea) ^\Tsa,)(^\Tea2), (6) 

as depicted in Fig. [5] (see Sec. II of SI). 

Figure [3] illustrates an intuitive description of the 
equivalence, which is mentioned in the first statement, 
between dynamics observed on a subsystem of a larger 
system. After defining the equivalence classes and the 
ME condition (in the second statement), we show that, 
once the ME condition is fulfilled, observing the natu- 
ral time evolution of SA, i.e., without active controls on 
them, is sufficient to specify the class (Theorem 1 in SI). 
The third statement is to ensure that even if the sup- 
port of psE in "Hb may move around, it is still possible 
to find a Hamiltonian Hse with which psE can be seen 
as a stationary state on SEA (Theorem 2 in SI). Then, 
this justifies the use of Eq. ([7]) below, or the mirroring 
effect of maximally entangled states, for our Hamiltonian 
tomography. 
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FIG. 3: An illustration of the physical situations that lead 
to identical observable dynamics. In this figure, the experi- 
menter measures the dynamics a small magnet (spin) shows 
by varying local control parameters. There could be multi- 
ple possibilities of external elements, i.e., environment, which 
would give rise to the same dynamics of the magnet, no mat- 
ter what he controls. For example, he cannot distinguish two 
situations; the whole laboratory may be in the magnetic field 
generated by wire that goes around the lab, or the external 
field of the same strength may be due to a strong perma- 
nent magnet. Similarly, in the case written in the main text, 
there could be many triples {(de, I'^'sea), Hse)} that lead to 
indistinguishable dynamics on SA. 



Now that we have two pairs of maximally entangled 
states, we move on to the Hamiltonian identification 
stage. One of the well-known properties of maximally en- 
tangled states is that an application of unitary operation 
on two subsystems at one side of the pairs is equivalent 
to that of its transpose on the other side. For I'i'sEA) in 
Eq. we can write 

UseI^^sea) = VaI^sba), (7) 

where Va = Ug^ acting on Ha- Therefore, the uni- 
tary evolution wc observe on the ancillary system A — 
A1A2 should contain information about the Hamiltonian 
Hse — i/thiUsEit). Naturally, however, simply looking 
at the state of A does not reveal any information on Use, 
but it turns out that identifying psA (t) by quantum state 
tomography suffices for our purpose. The origin of time 
t may be reset to the instance where the state-steering 
protocol is completed. 

Let us describe a specific method of estimating the 
Hamiltonian Hse after preparing the state I^^sea) in 
Eq. As all we can probe is systems A and S, we can 
only observe the time evolution of the reduced density 
operator p5A(i) = f''i^E\^SEA){'^SEA\. Our task is to find 
a Hamiltonian Hse that generates the time evolution of 
PSA{t), such that 

d 

i^PsA = [hJe ® Is, PS a], (8) 
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where H'g^ acts on T-La, despite its notation. The state 
PSA as a function of time can be obtained by iterating 
quantum state tomography at different times and is, thus 
a sort of inverse problem of the usual calculation of Eq. 
([8]). Note that in order to correctly specify H'g^ on the 
side of A as an interaction between two subsystems, the 
basis we employ needs to be written as a product of Ai 
and which are entangled with S and E, respectively, 
at the instance when the state \^sea) in Eq. ([6]) was 
realised. 

Since the system A is driven by Hj^, which is finite- 
dimensional, PS Ait) can be Fourier transformed such that 
each component pa evolves with a frequency 9a {a — 
0,1,..., i), where we set 9o = and L can be as large 
as dsdsidsd-E — l)/2. Namely, 

PsAit) = Po + E (P-^'"'' + Pi^"'^') ' (9) 

where all 6^ and Pa can be obtained from the data of 
PsA{t)- Substituting Eq. ^ into Eq. (|S|), we have 

[ffJij,Po] = (10) 

[H'sE^Pa] = daPa (H) 

[hIe.pI] = -OcpI (12) 

Now we choose a set of orthogonal bases of hermitian 

operators on T-La ® Hs as {/ij}jf,i^ so that any hermi- 
tian operators on the same Hilbert space can be written 
as a linear combination of hj. The relevant hermitian 
operators in our context are po, pa + Pa, i{pa — Pa) and 
^SE -^S- If we write the last one as 

HsE<»Is ^^Kjhj, (13) 

j 

Hamiltonian tomography can be seen as a task of iden- 
tifying coefficients kj. Expressing the above hermitian 
operators in terms of hj and using Eqs. p^ - p^ . we 
arrive at a system of linear equations with respect to Kj . 
The solution of the resulting set of equations would be 
non-unique. However, once we obtain a possible set of 
{kj}, it suffices to describe all the dynamics we observe 
through the accessible systems S and A, if any. Hence, 
the task of Hamiltonian tomography is completed. For 
the full description, refer to Sec. HI of SI. 

The Hamiltonian Hse thereby estimated contains all 
the necessary information to characterise the observable 



dynamics, albeit unmodulable per sc. What we can con- 
trol actively is the system S. Thus, the dynamics of the 
entire system SE is governed by the Hamiltonian 

H{t)=HsE + J2fn(^)Hs"\ (14) 

n 

(n) 

where He; are independent Hamiltonians that act on S 
and can be modulated by fn{t)- As we have already iden- 
tified Hse, there is sufficient information to judge the 
controllability of the system SE under this Hamiltonian 
P^ . A theorem from the quantum control theory states 
that the set of realisable unitary operations on the system 
is generated by dynamical Lie algebra [iGl - fTsj . Dynam- 
ical Lie algebra can be computed by taking all possible 
(repeated) commutators of operators in Eq. p4|) . i.e., 
iHsE and {zifg"''}, and their real linear combinations. 

Therefore, our knowledge of Hse allows the control- 
lable system to encompass not only the principal system 
S but also (a part of) the environment E. That is, we 
are now able to exploit the dynamics inside E for use- 
ful quantum operations, such as quantum computing, by 
controlling a small system S only. This is the same sit- 
uation as in refs. |19l423l |. where only a small subsystem 
is accessed to control a large system. 

We have shown the possibility of probing a large 
surrounding quantum system (environment) through a 
small principal system despite the lack of direct acces- 
sibility, provided the environment is effectively finite- 
dimensional. By probing, we mean fully identifying 
the Hamiltonian and utilising it as a useful resource for 
quantum control, e.g., quantum computation. This goes 
strongly against our intuitive supposition that environ- 
ment is usually a villain with respect to the protection of 
quantum coherence and that the Hamiltonian Hse is not 
estimable. Although our method might not be realisable 
in the lab straightaway, it opens up a path to the novel 
exploitation of the environment for quantum engineering. 
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Supplementary information: Probing an untouchable environment as a resource for 

quantum computing 



I. PROBLEM SETTINGS AND EQUIVALENCE 
OF DYNAMICAL BEHAVIOURS 

A. Problem settings 

We consider a joint system consisting of three parts: 
the principal system S, its surrounding system (environ- 
ment) E, and an anciUary system A, whose Hilbert spaces 
are denoted as Hs, "He, and Ha, respectively. The en- 
tire system SEA on Tis ®'He® T^a is a non-dissipative 
closed system. The principal system S interacts with 
its environment E via Hamiltonian Hse, while E does 
not interact directly with the ancillary system A. Thus, 
the Hamiltonian of the joint system can be written as 
Hse (Xi I A- We are not allowed to access the environ- 
mental system E directly, which means that no part of 
E can be a subject of direct control or measurement. 
Meanwhile, we are able to perform any quantum opera- 
tions and measurements on the joint system SA instan- 
taneously. 

A key assumption we make is that the environmen- 
tal system E is finite-dimensional, i.e., (Ie ■— dimi? < 
-l-oo. Those systems under our control, A and S are 
also finite-dimensional, and naturally dA ■= drmT-lA and 
ds '■= dim 7^5 are always known. We do not assume any 
prior knowledge ofdE, the interaction Hamiltonian Hse, 
and the state on Hsea ■— T^s ^ "He ® 'Ha- For most of 
the discussion in this supplementary information (SI), 
dA refers to the dimension of A after the state-steering 
protocol (see the main text and Sec IIVI in SI) has been 
completed, unless stated otherwise. 

Under these settings, our goal is to obtain as much 
information as possible about E and the interaction be- 
tween S and E, namely, Hei Hse, and the state on 
Use A- A central tool for the information acquisition is 
quantum state tomography of the joint system 5*^ to 
determine psA as a function of time. Hence, we assume 
that it is possible to prepare an identical, but not neces- 
sarily known, initial state \^sea) on the whole system 
'Hse A as many times as necessary. Pragmatically, such 
a state initialisation can be achieved by waiting for the 
equilibration of the state, which is caused by the interac- 
tion with a larger environmental system that surrounds 
i5 0, Q . The timescale for such an equilibration is much 
longer than the one within which the entire system of 
SEA can be considered closed. The initial time tg is 
defined as the time when the state initialisation is com- 
pleted. 

Since we have prior information about dA and ds, our 
system dynamics can be characterised by dE, a state on 
'Hse A at the time t^, and the Hamiltonian Hse- Note 
that at to, we can always assume the initial state on 
'Hse A is pure. This is because when the state on 'Hse A 



is mixed, we can append an extra Hilbert space T-Le' to 
the system so that the state on 'Hse A ®'He' is pure 
[1, Then, we simply redefine 'He^'He' as He, and 
Hse®Ie' as a new Hamiltonian Hse to restart the whole 
discussion. Therefore, we need a set of three elements 
{dE, \^ sea) , Hse), which we shall call a triple, to char- 
acterise the behaviour of our system under the effect of 
environment between the times and too {to < too)- Al- 
though we choose too = +oo, which is theoretically natu- 
ral, since we practically perform an experiment within a 
finite time length, the consideration of a finite too is also 
useful as we will see later. 



B. The equivalence due to indistinguishable 
dynamics 

As mentioned above, our primary goal is to identify 
the triple {dE,\'^SEA),HsE)- Yet, it appears to be a 
formidable, if not impossible, task to completely specify 
the triple when our access is limited to only S and A. 
What if there are more than one possible triple? In the 
context of system identification, the problem of indistin- 
guishable system models for a given experimental data 
set has been studied in the classical setting for a long 
time and also recently discussed in quantum scenario, 
in which the entire system is known to be controllable 

0. Similarly to these cases, for our task of identifying 
the system E, it turns out that even if there were multi- 
ple possibilities of triples that give rise to the same dy- 
namical behaviour on SA, the differences between them 
would not lead to distinct outcomes of quantum control 
^7] of E through S. In other words, if there were two in- 
distinguishable environmental systems characterised by 
{dE,\'^SEA),HsE) and {ds,]^ sea), Hse) in the time 
period [to,too), the results of any quantum computation 
that utilises E as a resource would be independent of 
whether the true environment was either of them. There- 
fore, for the acquisition of information on the environ- 
mental system T-Ie towards the exploitation of as a 
(partial) resource for quantum computing, it suffices to 
determine the equivalence class on the set of all triples 
through all possible sequences of quantum operations on 
'Hs^Ha- 

Before giving a rigorous definition of the equivalence 
class of triples, let us first specify all operations we can 
apply on the system. First, we do not consider opera- 
tions that are applied continuously in time. Thus, what 
we consider to be applicable is a sequence of instanta- 
neous quantum operations [1, 0, {r^j^^i at time 
ti, where n < +oo and ti < tj for all i < j. Sec- 
ond, a quantum operation Ti can be non-deterministic, 

1. e., trace non- increasing, because we can always post- 
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select the measurement results. Third, we are allowed 
to append and remove finite-dimensional ancillary sys- 
tems, which means that Ti is a quantum operation on 
S Hs) to <B {Ha, ® "Hs), where 03 {H) is a lin- 

ear space of all (bounded) linear operators on "H, HAi^i 7^ 
TiAi in general, dimH^i < +00, and Haq ■= Ha- 

Hence, our definition of the equivalence class is as fol- 
lows: 

Definition 1. (The equivalence between triples) A 

triple (dETl"^ sea) T Hse) ■5'^*'^ he equivalent to an- 
other triple, (dE.l'^sEA), Hse), in [to, too), if they sat- 
isfy 

Tte (jl (T.^Ie) o (Ia®U^%) PH'^sea))^ 

^Tte (^^®^e) o (lA^uj^l) PH^sea))^ (SI) 

for all n Cz N and all sequences of completely positive 
trace non-increasing maps {F,}"^]^ 0, 0, |j, 0/. Each Ti 
is a map from *B (Hyii_i ^'Hs) to %{T-LAi ^T-ls) and 
performed at time ti{i G {1, 2, n}). We will denote 
the equivalence between triples as {dE,\^ sea) , Hse) = 
[dE, \^sea),Hse)- 



In Eq. (jSlI) . I A and Te are identity (super)operators 
on T-La and He, respectively. Also, P(|5')) stands for 
P(|*)) := is given as 

^SsiP) ■= "^xp {~iHsE{ti - ti^i)) pexp{iHsE{ti - ^i-i)) 

a.nd Ug^{p) is defined similarly with instead of Hse- 
We can easily see that the relation "=" is reflective, sym- 
metric, and transitive. Thus, it is an equivalence relation 
in the mathematical sense [10[ , and a set of all the triples 
can be decomposed into equivalence classes accordingly. 

Our definition of equivalence here differs from the one 
in in that ours includes the possibility of appending 
an arbitrarily large ancillary system. Also, in Def. [1] 
above, the controllability of the system of interest is not 
assumed. Such a consideration is important especially 
when we can utilise the joint system SE as a part of a 
larger quantum network, rather than an isolated quan- 
tum computer. 

Let us slightly simplify the definition of the above 
equivalence relation for the following discussion. Here, 
we define a Hilbert space A' that includes all Ai as its 
subspace. 

Lemma 1. A triple {dE,\'^SEA), Hse) is equivalent to 
anther one {dE,\'^ sea) , Hse) in [^0,^00)7 */ and only if 
they satisfy 



SEA) 



n';1-u^%]\^sea) 



(S2) 



for all n € N, all sequences of real numbers {iijf^j^ with 
ti < tj for i < j , all finite dimensional Hilbert spaces 
Ha' which includes Ha o-s a subspace (Ha C Ha' ), o,nd 
all sequences of operators {N^sA'^^=i '"^ ^-4' ® ^S- 

In Eq. dSH), ■= exp[-iiJsB(t, -t^-i)], and Ie 

and Ia' are omitted. We shall not write identity oper- 
ators explicitly throughout the paper when there is no 
risk of confusion. 
(Proof) 

In order to prove the "only if part, since Ti{p) := 
Ng\,pNg^, is a trace non- increasing CP map, we sim- 
ply need to define Ha' as Ha' '■— HAi- For the "z/" part, 
we define Ha' '■= ®"=o^^i- Then, the linearity of the 
partial trace Tr e guarantees the statement of the lemma. 



II. THE MAXIMAL ENTANGLEMENT 
CONDITION AND THE EQUIVALENCE 
BETWEEN TRIPLES 

We now introduce a condition that is of crucial im- 
portance for our analysis. We shall refer to it as the 
maximal entanglement (ME) condition. It will be shown 
that when our system satisfies this condition, the natural 
time evolution of the system SA completely determines 
the equivalence class (the subsection III A[) . By natural 
time evolution, we mean the evolution of the system with- 
out active operations on it, namely, the evolution that is 
driven only by the system Hamiltonian, which is Hs e in 
our case. Then, we prove that when the system satisfies 
the ME condition, there exists an equivalent triple, in 
which E which is maximally entangled [1, [ll|, [13] with (a 
subset of) A (the subsection III B I) . We note that although 
our main focus is on a finite-dimensional in this pa- 
per, some of the theorems and lemmas in this section are 
valid even for infinite-dimensional He, b^s long as dA and 
ds are finite. 



A. The mELximal entanglement condition 

The condition for the equivalence of triples still ap- 
pears quite complicated even in the form of Lemma [T] 
In this subsection, we prove that the condition for the 
equivalence reduces to merely the indistinguishability of 
the natural time evolution of the system SA, when the 
condition defined in the following is satisfied: 

Definition 2. (The maximal entanglement condi- 
tion) The maximal entanglement (ME) condition is said 
to be satisfied by a triple {dE,\'^SEA),HsE) in [^Oi^oo); 
if for all t G [^Oi^oo), there exist Hilbert spaces HAi{t) 
and Ha2 (t) such that Ha — Hai (t) <E) Ha2 (t) , dim Hs — 
dimHAi{t), and 

TteP il^sEAit))) = P {\TsA, m (E> PA, (t), (S3) 
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where \TsAi(t)) is a maximally entangled state onHs ® 
HAiit) and a state PA2ii) ^ projector onto HAii^) (up 
to a proportionality constant). 

Throughout this paper, \T xy) denotes a state that 
is maximahy entangled fully on the space T-Lx ® 
Hy specified by the subscripts. That is, \Txy) = 
^x^^^Sffi K^)Nv), where {\ix)} and {liy)} are or- 
thonormal bases of Hx and "Hy, and dx = dimHx{— 
rfy). 

In the above definition, \'^sEA{t)) is the entire state 
at time t, i.e., 

I^-sbaW) - exp(-iiJsB(i - to))\-^sEA)- (S4) 

Here, we note that the Hilbert spaces Hai {t) and 'Ha2 {t) 
may vary inside T-La as the time evolution of SE (due to 
Hse) would be reflected in Ai and A2 through entangle- 
ment. 

Equation (jS3p implies the existence of a pure state 
l^EAiit)) on He ®'HA2{i) such that 

\^SEA{t)) = \TsA^{t)) ® \^EAAi)). (S5) 

where |$_Eyi2 (^)) is a maximally entangled state (MES) on 
a subspace oi^E® H-a^ (t) ■ Yet the rank of the reduced 
density matrix pE may be smaller than the dimension of 
the system E: rankp^; < d^. It turns out, however, it is 
possible to choose a triple {dE, \^sea),Hse) (equivalent 
to (d^;, I^Psba), -ffss)) so that \^EA2{t)) can be expressed 

as \TEA2} = d'E^^^ J2i=i Na2)N-b), i-e., a state that is not 
only maximally entangled but also satisfies rankpg = dE 
(see Theorem [5] and Corollary [3]) . When the ME condi- 
tion is found to be satisfied, we can take it for granted 
that the whole E is maximally entangled with A2 C A, 
despite the inaccessibility of E. 

An observation is that the fulfillment of the ME con- 
dition can be tested through tomography of the state on 
SA only. This fact leads to a lemma: 

Lemma 2. Suppose {dET\'^SEA),HsE) satisfies the 
ME condition in [ta,too), and {dE,\^ sea) , Hse) = 
{dEA'^SEA),HsE) in[to,toc)- Then, {(IeA'^ sea) , Hse) 
also satisfies the ME condition in [to, too). 

(Proof) Due to the definition of the equivalence, if 

{dE,\'^SEA),HsE) = {dE,\^SEA),HsE) in [to, too), we 
have, for all t e [to, too), 

Ti-eP il'i'SEAit))) = TveP (l^SEAit))) , (S6) 

where \^sEA{t)) is given in Eq. (|S4|) and |^s_ea(0) is 
defined similarly. Then, since the fulfillment of the ME 
condition only depends on the reduced density opera- 
tor on Hs ® T-La during [to, too), Eq. ((S6|) guarantees the 
statement of the lemma. I 
Lemma [2] implies that the ME condition can also be con- 
sidered as a property of equivalence class of triples. 

Let us now present a theorem, which claims that if two 
triples give rise to an identical (natural) time evolution 




FIG. 4: Schematic illustration of the relation in Eq. (|S9|l . 
System A' includes Ai and A2 as its subsystems. No matter 
what operation Nsa' is performed on SA', the same effect 
can be realised by an operation J-p^^^ {Ngj^i ) applied solely on 
A'. 

on SA and if one of them satisfies the ME condition, then 
those two triples are equivalent. 

Theorem 1. Suppose {dE,\'^SEA),HsE) satisfies the 
ME condition in [to, too)- Then, {dE,\'^SEA),HsE) = 
(d-E, l'^ sea), Hse) in [to, too), if and only if they satisfy 

TteP (|*sBA(i))) - Tr^P {[^sEA{t))^ (S7) 

for all t e [to, too)- 

Thus, when the ME condition is satisfied, an identical 
time evolution on SA is sufficient to certify the equiva- 
lence, that is, we do not need to consider all possible se- 
quences of quantum operations {Ti}. Towards the proof 
of the theorem, we show two lemmas. 

Lemma 3. Suppose a state ["^sea) on the Hilbert space 
Hs ^He ®y-Ai ® 'Ha2, where ds — dim'Hs = dim'HAi, 
and ["^sea) can be written as 

[■^SEa) = |T5aJ ® [^EA2), (S8) 

where [TsAi) is a maximally entangled state, and 
TteP {[^eA2)) *s projector up to a proportionality con- 
stant. Then, for all finite- dimensional Hilbert spaces Ha' 
that include Hai ® 'HA2 as a subspace, and for all oper- 
ators NsA' on Hs ® Ha' , 

1. su^v^iA'P {NsA'['i> sea)) C suppTr^Pd^-sisA)) 

2. For a given state psA '-= '^'^eP {\'^ sea)), there 
exists a linear map J^psA from ^{'Hs®'Ha') to 
58 {Ha') such that 

J" PS A {NsA') ® Ise['^sea) = NsA' ® Ie['^sea), (S9) 

for allNsA' e^iHs^HA')- 

Here, we note that the same Tp^A satisfies Eq. (jS9p 
for all states ['9 sea) having the common reduced den- 
sity matrix psA = Tr^P (I^'sea))- Figure |4] depicts the 
equivalence relation of Eq. ([59]) . 

(Proof) 

The proof of the first statement proceeds as follows. We 
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write \TsAi) and \^eA2) as 

1 



1$ 



1 



EA, 



ds 

E 

r 



where {le.)}:'^,, {I/,)};;,, and {b)};:r'^ are 

orthonormal bases of 'Hs^T^Ej'Haa, and Ha2j respec- 
tively, and r is the Schmidt rank [3| of \(^eA2)- Without 
loss of generality, we see 



, dim 'H / 



NsA' ®Ie\'^sea) = 



R 

E 

a=l 



|?7a)sB, (SIO) 



where { | ^„ ) j)^™ and { 1 77^ ) } 



dsdB 



of T-La' and T-Ls 



a=l 

"He, respectively. 



are orthonormal bases 



Also, R and 



{Xa}c,=i are the Schmidt rank and the Schmidt coef- 
ficients [3] of Nsa' (8> IeI'^sea)- Naturally, equation 
TiSA'PiNsA'l'bsEA)) = TtsaPH'S'sea)) holds, and by 
substituting the above expressions into each side we have 



R ^ 

^"^S^KMiValsE = - 



(Sll) 



which means suppTr5|77Q)(77Q| C span for all 
a. Then, \ria) should be in the space spanned by 



{\(^^)\fj)}t=ij=i' which in turn implies span {| 770)}"^ 1 C 
span{|ei)|/j)}^^^ ^^j^, hence statement 1. 

We now move on to the proof the second state- 
ment. Suppose {|ga)}^='i is an orthonormal basis of 
Ha': where dj^i = dimH/i'. There is a natural lin- 
ear isomorphism, /, from Ha' ® suppTr^jP H'^sea)) to 
^ {suppTr aP il"^ sea)) ,'Ha') (*B(-H,/C) is a space of ah 
linear operators from H to /C). / is defined by the corre- 
spondence of their bases: 



R 



f ■■ l5o) 



l/,)^l5a)(e.|®(/j| (Va,*,j-). (S12) 



Then, we define a linear map J^g^ as 

J'psANsA') := Vd^f{NsA'\^SEA))V, 



(S13) 



where Nsa' G S {Hs (8) T-La'), and y is a partial isometry 
defined as 

ds r 

y--=J2T.\^^^s\fj)EmAAj\A2. 

j=i 

We note that J-p^A is well defined thanks to statement 
1. Then, it is straightforward to see that J^ps^ satisfies 
Eq. dMl)- 

It is possible to show that the linear map J-pg^ only 
depends on a reduced density of matrix psA of the state 
1 5*); more specifically, its effect does not depend on the 



Schmidt bases {|ei)}fji and Although / and 

V do depend on {|e,)}ffi and their depen- 

dence cancels out in Eq. (jS13|) . As a result, any specific 
choice of these bases does not affect the action of J^ps^ ■ 
We can also confirm this property of Tps^ from the fact 
that equation J^p^^ (Nsa') C/bI^sba) = Nsa'Ue]-^ sea) 
holds for all unitary operators Ue on He- Therefore, 
statement 2 holds. ■ 

Lemma 4. Suppose {(Ie, \^sea)tHse) satisfies the ME 
condition in [toi^oo)- Then, for all instances {ti\^^l, 
where to < ti < tj < too for i < j , all finite dimensional 
Hilhert spaces T-La' 3 T-La, cind all sets of linear operators 

|-^SA'} . ""^ ^A' , the following equation holds: 



n 

n H^c^^) i'^sea) 




J=2 



(j) X- 
SA'-^PSAit 



SEA 



(tn)), 



^SE 



(SI4) 



where NsA'^p is a linear map on *8 (T-Ls ® "Ha') defined 
as NsA'J'piMsA) = NsA' {Tp{MsA') ® Is) for Msa' e 

(Proof) 

By repeatedly applying Lemma [31 we have: 



n 

n {^sMe 

3 = 1 

n 

n {^siui^i 

n 

n {^siu^^E 

n 

n {^sWPe 

n 

n i^sMi 



1^ 



SEA/ 



N^ll^SEAih)) 

{n^sI) I'^'SEAih)) 
^SA'^PSA(ti) (iVi'i) \^SEAit2)) 



PSA{t 



3=3 



t.) {N^sl'^psAiu) {n^sI)) 




iSE 



|*SBA(i«))- 

Thus, we derived Eq. (|ST4)) ■ 

We are now ready to prove Theorem [TJ 
(Proof of Theorem [l]) 

The "only if part is trivial from the definition. We 
therefore prove the "i/' part of the statement. Sup- 
pose (d_E, I^'sba), -ffsfi) satisfies the ME condition and 
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idEA^SEA),HsE) and idE,\^SEA),HsE) satisfy Eq. 
((57|) . Then, for all instances {ti}^^^, where to < ti < 
tj < too for i < j, all finite-dimensional Hilbert spaces 

Ha' 3 Ha and all sets of linear operators l^/V^^/ 1 on 

® Ha', we derive the following equations: 



(I'i'SEAitn))] 

71 



TteP 





FIG. 5: A situation in which the ME condition may be satis- 
fied, but the support of the state \^eA2) in E moves around 
in a larger space. Theorem [5] rules out such a structure once 
the ME condition is fulfilled. 



complex numbers z. Together with Theorem [TJ we reach 
the statement of the corollary. ■ 



i=2 



i=2 



Tr EP{\^SEA{tn)) 



--TveP n {n^a'UsI) I^sea) 



t 



where Eq. (|S14p of Lemma |4] is used in the first and third 
equalities. The above equation is nothing but Eq. ((S2|) 
in LemmalU and thus the "if part of the theorem holds. 
■ 

Theorem [T] leads to the following corollary. 

Corollary 1. Suppose {(IetI'^sea), Hse) satisfies the 
ME condition in [to, toe), o-nd ds < +oo. Then, 

idEA^SEA),HsE) = idE,\^SEA),HsE) in [to, too), 

and dE < +oo imply that {dE,['^SEA{io)),HsE) = 
[dE, \^SEA{to)),HsE) in any time interval [iojoo)- 

Roughly speaking, what this claims is as follows. Sup- 
pose there are two triples, the one with E and the other 
with E. If the ME condition is fulfilled by one of them 
and both give rise to the identical natural time evolution 
on 5*^ for a finite duration, regardless of its length, then 
both triples are equivalent; that is, any active quantum 
operations on SA cannot reveal the difference between 
them. 

(Proof) 

Since both dE and dE are finite, both T^tieP {\^ SEA{t))) 
and TteP [[^ SEA{t))^ are analytical functions and 
coincide on [^oj^oo)- Thus, by analytical continuation, 
1:veP[[^SEa{z))) = TrBP(|«'sB^(0))) holds for ah 



B. The ME condition implies maximal 
entanglement between A and SE 

In the previous subsection, we have suggested that 
when the triple [dE, ["^ sea) ,Hse) satisfies the ME con- 
dition in [to, too), there exists an equivalent triple in 
which I '^EA2 (t) ) can be chosen as a maximally entangled 
state on the entire space oi'H.E^'HA- Now we shall prove 
it as a theorem in this subsection. 

Theorem 2. Suppose {dE,[^ sea) , Hse) satisfies the 
ME condition in [to, too), o-nd dE < +oo. Then, 
there exists a triple {dE,['^ sea) , Hse) such that 
{dE,['^SEA),HsE) = {dE,['^SEA),HsE) in [to, too), and 
["^SEa) is a maximally entangled state on the full space 
of SEA with respect to the partition between SE and A; 

that is, TtaP {[^sea)) = djd^Is «) Ie- 

Let us sketch the idea of its proof: For 
{dE,['isEA),HsE) satisfying the ME condition, 
['^SEa) is a maximally entangled state on the full 
space of S,E, and A if and only if dE — rank(0£;(t), 
where PS (t) := Trg^^^ (e"*^^^^*-*") |*s_ea)) • If the 
rank of pE{t) is smaller than dE, suppp£;(i) may move 
around within E as time proceeds (see Fig. [S] for an 
intuitive illustration). In such a case, the property 
of maximally entangled pairs, Eq. (7) in the main 
text, cannot be used for our Hamiltonian identifica- 
tion purpose. Thus, in order to prove the theorem, 
we need to show the existence of Hamiltonian H'g^ 
such that {dE,['^sEA),HsE) = {dE,['i> sea),H'si^) 
and [H'g^, psE{to)]=^- With such a H'se, 

p'j^it) := Tr5AP(e-*^si=;(*-*o)|^^^^)^ be time- 

independent. 

To this end, we restrict the space He to suppp'^{t), or 
equivalent ly, we define dE, ['^sea), and Hse as dE '■= 
lankpEito), [^sea) ■= I^sba), and Hse ■= {IsA <^ 
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PE{to))H'g^{IsA<^ PE{to)), respectively, where PE{t) := 
rank/9^ (t) • p'^{t) is a projector onto the support oi pe- 
Then, the triple (d^, \^ sea)tHse) satisfies the desired 
conditions: {dE,\^ sea) ■, Hse) = {dEA"^ sea),Hse) in 
[to, too), and \^sea) is a maximally entangled state on 
the full space of SEA. Therefore, the proof of the theo- 
rem can be reduced to constructing Hamiltonian H'g^ 
so that {dE,\'^SEA),HsE) = (de, I^sb^i), ^^ss) and 
[Hg^, psE{to)]=0- In fact, we shall see that the ME con- 
dition enables us to do so. 

In order to construct the H'g^, we need more detailed 
description of Hse- In general, a Hamiltonian Hse can 
be decomposed as 



Hse = /s ® Hid 



where {IsjUWcjoLi^ 



4-1 



Hn 



(S15) 



forms an orthogonal basis of a 
real space of all Hermitian operators on Hs- When 
{dE, I'i^ sea) , Hse) satisfies the ME condition, the state 
on SE is proportional to Is <E} PE{t), where PE{t) is a 
projector-valued function on E. Then, the Schroedinger 
equation for \^ sEA{t)) becomes 



ij^Is®PE{t) 

--[HsE,Is®PE{t)] 



(S16) 



--Is ® [Hid,PE{t)] + ^ <T„ $5 [H^,PE(t)] , (S17) 



for alH e M. Thus, H^^ and H^ satisfy 



'iJ^Eit) = [Hid,PE{t)], 

dt 



Q=[H^,PE{t)] (Va), 



(S18) 
(S19) 



That is, only Hid is responsible for the time evolution of 
PE{t) = -PE(t)/rankp£;(to)- Therefore, a modification to 
-ffid may lead to desired Hamiltonian H'g^. Suppose that 
the desired modification is described as Hid Hid — H[^ 
with a Hermitian operator H[^. Such a H[^ is shown to 
exist and has properties as in the following lemma. 

Lemma 5. Suppose He is a finite- dimensional Hilbert 
space, and suppose a projector-valued function PE{t) on 
T-Lej and Hamiltonians Hid '^^'^{^ala^i onH.E satisfies 
Eqs. I57g|) and (ST^ for all t € It, where It is an open 
subset of M. Then 
satisfies the following equations 



there exists Hamiltonian H'^ that 



[Hid,PEit)] - [Hld^PEit)] 
[Hid,Hl^] = 0, 
[H^,H:^]=0 (Va). 



(teit. 



(S20) 
(S21) 
(S22) 



(Proof) 

We generate a sequence of triples | Gn, H^^^^ | 



(n) 

where H^^ are Hamiltonians, Qn are linear Lie algebras 
on T-Le, and G„ are Lie groups [l^, according to the 
following procedure. Suppose Ha have spectral decom- 
positions as Ha — J2j ^l"'^ Pj"^ ' where h'j"^ ^ /i^?"* for all 
j 7^ j' , and Pj"'' are all projectors. Then, gi is defined as 
a Lie algebra whose elements are (Vj, a), Gi is 

a compact Lie group generated by Qi, and H^^ is defined 



as 



^id - 



UHidUUp (U) , 



(S23) 



ueGi 



where p{U) is a Haar measure [T^ [TJ] on Gi. 
Starting from (^gi,Gi, H^^^ , we recurrently define 



0„, Gn,H- ?^ ] as follows: Suppose Hid 



-H. ?^ has a spec- 



tral decomposition as Hid — ^m'' = ^jQj 



where 



are projectors. Then, we let Qn+i and G„+i be a 
linear Lie algebra consisting of (^j) ^^^^ 

a compact Lie group generated by Qn+i- Similarly as 



above, H^^'^^^ is defined as 



rr("+l) . 

-"id 



UeG„ 



UH[2^uUp{U). 



(S24) 



Then, we can prove the following equations for all t d It 
(The time t will denote t £ It hereafter unless otherwise 
stated): 



[Ha,Hld^]=0, (Va,n) 

[h,-hI-\hI-+'^] = o, 

[Hid,PEit)] = [Hl'd\PEit)l (Vn). 



(Vn) 



(S25) 
(S26) 
(S27) 



Note that Eq. (|S37[) can be reexpressed as [Qj , Pe {t)] = 
by subtracting the RHS from the LHS. 

Let us first show the following two equations: 



and 



^] '-"id 



h '^id 



0, (Vj,a,n) 



(Vj,n), 



(S28) 



(S29) 



from which Eqs. (|S25p and (|S26|) are directly obtained. 



Since Us :== /b -P}"^ -f e^'^P)"-' is in G„ for all 5 € [0, 27r], 



iS nia) 



h[2^ satisfies 



H 



(n) 



id 



5G[0,27r] 



/-P 



P 



(a) 
(«)/ 



-j5p(a) 



2^ 



id 



p(a) jj(n) p(a) 
-"id ' 
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which imphes 



'id 



Thus, Eq. ([S28l) holds. Similarly, we can prove Eq. (fS29l) 
from the fact Us :^ Ie — Q'j^^ 

We now prove Eq. ((S27l) by induct ion. First, we con- 
sider the case n = 1. Equation (|S19[) implies that pj"-* 

f\PE{t)] 



satisfies 



e'^Q(") in G„+i. 



0. Then, by differentiating this 



equation and using Eq. (|S18I) . we obtain 







(S30) 



for all j and a. Then, because of the following identities. 



p^"\[i/id,PB(t)] 



0, 



p(ai) p(Q2) 
il ' 12 



,P 



p("3) 



p("2) p(Q3) 
12 ' is 



0, 



,p; 



("i) 



the last of which is the Jacobi identify, all operators writ- 
ten in the form of 



p(ai) p("2) 



("3) 



,p: 



(S31) 



commute with Psit) and [i?id, PE(t)]. Thus, since an 
arbitrary X G gi can be written as a linear combination 
of the terms in the form of Eq. (|S3ip . any X £ Qi satisfies 
[X,PE{t)] = and [X, [H,d,PE{t)]] = 0. 

Since Gi is a compact Lie group, for a given U £ Gi, 
there exists an operator X G gi satisfying U = expX. 
Thus, an arbitrary U Gi satisfies 



and 



[U,PE{t)]^0 



[C/,[ffid,P£(t)]] = 0. 



(S32) 



(S33) 



Hence, by using Eqs. (jS32p and (|S33p . we can derive the 
following equation for all U E Gi: 

[Hid, Psit)] = U [Hid, PEit)] C/t 
- [UHidU^,PE{t)] . 

By integrating the above equation on Gi with a Haar 
measure n{U), we obtain [Hid, PE{t)] = H[^\PE{t) ■ 

Second, following a discussion similar to the case of 
n = 1, we prove Eq. (|S27I) for n -f 1 assuming that Eq. 



((527| holds for all k < 



k < n. As in the case of Pj"\ this fact implies 



Qf\PE{t) 



for all j. 



Qf\[Hid,PE{t)] 



= 



(S34) 



for all j, k < n. Then, since with all A: < n and 

(iPj"^} generate G„+i, Eqs. ((S30)) and ((S34l) imply 

that an arbitrary X e 0n+i satisfies [X, PE(i)] = and 
[X, [Hid, Psit)]] ~ 0. Therefore, as in the case of n ~ 1, 



we deduce [Hid, PE{t)] = 



H 



(n+l) 



id 



PEit) 



Equations. ((S25)) . ([S26l) . and ((S27| can be used in 
the final step to prove the lemma. Since all 0„ are Lie 
sub-algebras oi su{diui'HE), and thus finite-dimensional, 
there exists an TV S N such that 0„ = Qn+i for all n > 
(note Qn C 0ri+i)- Therefore, we have 



[Ha , H.^ 
[Hid 



id J 



0, 



„{N) „(Nh 
^id '-"id J 



0, 



[Hid, PEit)] = [H[r,PE{t)]. 



(S35) 
(S36) 
(S37) 



Taking :— H^^^ proves the lemma. ■ 

Let us proceed to the proof of Theorem [2l 
(Proof of Theorem [2]) 

Suppose {dE, I^'ssa), -ffsfi) satisfies the ME condition in 
[to, too), i.e., PE{t) tTysaP {['^SEA{t)}) is a projector- 
valued function, where r := dimsuppTr5/iP (|^'5e^)). 
Therefore, Eqs. ((gT8)) and ((gT9)l hold with Hid and 
{HaY^^i defined in Eq. (jS15[) . Lemma [5] then tells that 
there exists a Hcrmitian operator H^^ on He that satis- 
fies Eqs. (IS20l)-(IS22ll for aU /* C M. We now define H'^^ 



H'qt^ :— HsE — Is ® Hi, 



(S38) 



Then, [H'sj^,Is®H!J and [H'^j^, Is ^ PE{t)] can be 
computed as 

[H'sE,Is ® HU] ^Is ® [Hid - HU, HU] 

=0 (S39) 



and 



[H'sE, Is ® PEit)] =Is® [Hid - HU, PEit)] 

= 0, (S40) 



where we h ave u sed Eqs. (|S2ip and (|S22p in Eq. (|S39p . 
and Eq. ((520)) in Eq. ((S40)) . Then, by defining 
\%EAit)) '■= exp i-iH'sEit - to)) I^-sba), we obtain 

TTEPQ'^ifsEAit))) 

=TrE exp i-iHsE{t - to)) P (I^-sea)) exp iiHsEit - to)) 
^TvE exp i-i ils «) F/j) it - to)) exp i-iH'g^it - to)) 

■Pi[-9sEA))eJip ilHsEit-to)) 

•exp ii ils ® HU) it - to)) 
=TiE exp i-iH'sEit - to)) P (|*5ea)) exp iiH'sEit - to)) 

^TTEPi[%EAit))), (S41) 



where we have used Eqs. (|S38p and (jS39p in the second 
and third equalities, respectively. Also, since the state 
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on SE can be written as 

TrAP{WsEAm 

■ {Is <E) PE{to)) exp{iH'sE{t - to)) 

^^-^s^^^bIM, (S42) 
rds 

the state TtsaP (I'^'sEAi't))) = 7-Pe(*o) does not depend 
on time t. Therefore, by defining a Hilbert space He, 
state |\&s_ea), and a Hamiltonian Hse on Us ^He as, 

-He := suppP£;(io), (S43) 
I^sea) ■■= l^EAito)) = I^sea), (S44) 
Hse := Is «> PE{to)H'sEls ® PE{to), (S45) 

the triple (^(IeiI'^sea), Hse^ satisfies the statement of 

the theorem, where dE dimHs = rankPE(io)- B 
An important corollary of Theorem [2] follows: 

Corollary 2. Suppose {dE,Y^ sea) , Hse) satisfies the 
ME condition in [ioj^oo); o,nd diva dE < +oo. Then, 
{cIeA^ SEA{t'o)),HsE) with |*5£;A(to)) being given by 
Eq. (S4\ ) also satisfies the ME condition in [^0,^^) for 
all t'^ and t'^ {t'^ <t'^). 

(Proof) 

Theorem [2] guarantees that there exists 
{dEA^SEA)-,HsE) such that {dE,]"^ sea) ,^Hse) = 
{dE,\^SEA),IIsE) in [to, too) and \^sea) is 
a MES with respect to the partition A\SE. 
Then, for all tg and t'^ satisfying t'^ < t'^, 

{dE,\^SEA{t'o)),HsE) = {dEA'^SEA{t'o)),HsE) in 
[^0!^'oo)i ^-s stated in Corollary [T] On the other hand, 
since \^sea) is a MES, idE,\^SEA{to)),HsE) satisfies 
the ME condition in any time intervals [Iq, t'^) (tg < t'^). 
Recalling the definition of the triple equivalence (and 
also the "if part of Theorem [T]), this leads to the 
statement of the corollary. ■ 
This corollary tells that if our system satisfies the ME 
condition for a non-zero time period, no matter how 
short it is, it will always satisfy the ME condition from 
then on. 



III. TOMOGRAPHY OF THE ENVIRONMENT 
UNDER THE ME CONDITION 

In this section, we show that, when the entire sys- 
tem SEA satisfies the ME condition, it is possible to 
specify the equivalence class by performing tomography 
on joint system 5*^. More precisely, when our triple 
{dE, I^'sba), Hse) satisfies the ME condition in the time 
period [to, too), we can determine dE, \^sea) and Hse 
that satisfy [dEA"^ sea) , Hse) = (dE,\^SEA), Hse) in 



[to, too)- As we have explained in Sec. U this informa- 
tion is sufficient to control system E and exploit it as a 
resource for quantum computation. 

Let us start with another corollary of Theorem [51 

Corollary 3. Suppose a triple {dEjl^^ sea) , Hse) sat- 
isfies the ME condition in [to, too)- Then, there 
exist Hilbert spaces Hai o-rid Ha^ o,iT-d a triple 
(^dE,\^ sea),Hse^ , such that Ha = Hax ®'Ha2, o.'^d 
I'^SEa) can be written as 

\-^sea) ^\TsA^)(^\rEA,), (S46) 

where \TsAi) G Hs <E) Hai and \TEA2) € 'H-e ^ 
'Ha2 are maximally entangled states, and dE — 
dimsuppTr^2^ (I'f -EA2))- Moreover, we can choose a ba- 
sis set freely for both A2 and E such that \T EA2) can be 
expressed as 

\TEA2)^^Y.\iA2)\^E). (S47) 

VdE 1=1 

The last part of the corollary states that there is free- 
dom in the choice of a basis for A2 and E due to the 
equivalence of observable dynamics induced by Hamilto- 
nians that may be equivalent up to a local unitary Ue- 

(Proof) 

The first part of the statement is simply a rephrase 
of Theorem [21 Thus, we only prove the second 
part, which says that IT^;^^) can always be writ- 
ten in the form of Eq. (IS47[) even if we choose 
arbitrary bases. Suppose that idEjl"^ sea) , Hse) 
satisfies the ME condition in [to, too) and \^sea) 
and Hse are given as \^sea) ■= {I <E) C/£;)|*sba) 
and Hse ■= {Is Ue)Hse{Is ^ U^), respec- 
tively, where Ue is an arbitrary unitary opera- 
tor on He- Then, \^sEA{t)) ^ := e-'"^^'\^ sea) 
and \'^SEA{t)) := e-*^^^*|*sBA) satisfy 

TTEP{\'^SEA{t))) - TrE{\^SEA{t))) (Vt G [to, too), 

i.e., {dE,\^sEA),HsE) = (^b, I^sba), -ffsfi) - Hence, 

the state \'^sea) has a unitary freedom on E, i.e., the 
freedom in choosing a basis to express it as Eq. (jS47p . 
■ 

Therefore, when our system satisfies the ME condi- 
tion, we can determine the decomposition T-La — Hai ® 
Ha2, the state \TsAi), and the dimension dE satis- 
fying dE = dimsuppTrs_B^jP(|5's£;A)) by performing 
(joint) state tomography on SA. Moreover, by redefining 
Ha ■— suppTr seP {\'^ sea)), we can assume dimH^^ = 
dim^B = dE and dim "Ha — dsdE- The above corollary 
also implies that we can always assume system EA2 is in 
the dE X d_B-dimensional standard maximally entangled 
state (of the form of Eq. ((S47)) ). 

The remaining task is now to determine the interaction 
Hamiltonian iJsB (finally!). Theorems [T] and [5] allow us 
to state the following: 
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Corollary 4. Suppose a triple {cLe, \^sea),Hse) satis- 
fies the ME condition in [to, too)- If o, Hermitian matrix 
HsE on Hs ®'He satisfies 



Is 



(S48) 



for all t in a neighbourhood of to, then 

{dE,\^SEA),IIsE) = {dE,\^SEA),HsE^ in [to, too)- In 

Eq. (S4§, psA{t) := Tr£P(e-^^s«(*-*°)|^'sBA)): the 
transposition T is taken with respect to the Schmidt basis 
of ['^sea) and Hg^ is an operator on Ha- 

Thanks to Corollary 3, we can take it for granted that 
['^SEa) has the form of Eq. (|S46I) . Corollary H implies 
that, for a given psA{t), which can be specified by state 
tomography, all Hermitian matrices HgE satisfying Eq. 
(jS48|l can be adopted as the interaction Hamiltonian be- 
tween systems S and E. 

Let us now describe how to find such a matrix Hse 
from the observed data of psA{t)- Since the time evo- 
lution of psA(t) is induced by the (finite-dimensional) 
matrix Hsei there exist a set of real numbers {0a}a^i 
and a set of linear operators {pa}a=o '^^ ® such 
that psAif) can be written as 



PSA 



(<) = po + Y. (e''°^'~*°Va + e-^'-'^'~"^^pl) , (S49) 



ct=l 



where L is at most dsd-E (^dsdE — 1^ /2 and po is Hen 
tian. Setting t = to, we have 



PSA{to) = Po + Y{pa+ pi) 

a=l 



Further, differentiating Eq. (IS49P n times leads to 



if and only if psA{t) € *B [T-Ls ® Ha) given by Eq. f5^gp 
satisfies 



[H, PS Ait)] 



(S53) 



In Eqs. (|S52p and (IS53[) . Is is omitted for simplicity; 
that is, H in these equations means H ^ Is- 
(Proof) 

Simply taking the time derivative of Eq. (jS49|) and then 
using Eq. (|S52[) lead to Eq. (jS53p to prove the only if 
part. The i/ part can be shown by substituting Eq. (|S49p 
into Eq. §53\l- 



1j 

[H,po\ + Y,{i\H,Pc.]-Oo,Po:)e 
+ {[H,pi]+e^pl)e-^'-'] 



= 0. 



(S54) 



Equation (|S52p follows if we apply 
liniT^+co T /o^'^^e*''"'* to Eq. (IS54l) . ■ 

Lemma assures that any Hermitian matrix H satis- 
fying Eq. (jS52p can be identified as an Hse (taking its 
transpose H'^). The existence of such a matrix H, i.e., 
a Hermitian matrix that satisfies Eq. (jS52p . is guaran- 
teed by Corollary m Therefore, our task now is to find a 
specific form of H, given the information on psAit)- 

To this end, we fix an orthonormal basis of all Her- 
mitian operators {Aj}'^.^^'^ c ^ {Hs 'E) Ha), each of 

which has the form Aj = Bj (g) Is, where {Bj}^^-^'^ C 
55 (Ha) (Recall that our target Hamiltonian has the form 
Ha (8) Is and pAiS is maximally entangled.). We also let 
(S50) {ejfe/},7j; {j,k,l £ {1, 2, d|c?|;}) denote a set of struc- 



ture constants of the Lie algebra generated by , 
i.e., 



-psAit) =^{(z0„)'V„ + (-z0„)"pt}. (S51) 

t—to 



a=l 



Therefore, we can determine {0a}a=i and {pa}a=o from 
at most d^-th order derivative of psA{t) ai t = to, which 
can be obtained experimentally in principle. 

The information on {Oa\^^i and {pa\^^o aUows us to 
reconstruct Hse so that it satisfies Eq. (IS48p . Here is a 
lemma that shows a property the desired matrix should 
have: 

Lemma 6. A set of real numbers {Qa}]^^i, a set of lin- 
ear operators {/Oq}^^q C *B {Hs ® Ha), and a Hermitian 
matrix iJ e 05 {Ha) satisfy the following equations: 



[H, Po] = 
[H,pc.] = ~ 



9^Pa, (l<Va<L,; 



(S52) 



4dl 

[iA.,iAk] = e^kiiA, 



1=1 



(S55) 



Since the basis set is a basis of all Hermitian 

operators on Hs ^ Ha, we can uniquely expand the fol- 
lowing operators in terms of {^j}: 



Po = uiA 



4d| 



.7 = 1 



i {Pa - pi) = 



(S56) 
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where , ^Vj"''^ | , and | are all real 

constants for all a. Similarly, an arbitrary Hcrmitian 
operator H on Ha can be written as 



(S57) 



with real constants {ft-jj^f^'J^ ■ Then, a necessary and suf- 
ficient condition for H , {9a}^^i, and {pa\^^Q to satisfy 
Eq. (|S52I) is that {hj}^^-^'^ satisfies the system of linear 
equations with 1 <l < dlgCp^ and \ < a < L: 



2 j2 / ,4 J2 




^ ^ ejfc/Wfe I /i," = 0, 



(S58) 



The system of linear equations (IS58P has at least one 
solution, as for experimentally obtained {da}a=i ^'^'^ 
{Pq}q=o there must exist H that satisfies Eq. (|S52|) . In 
fact, the solution to the above equations, Eqs. (|S58|) . is 
not unique, since there are dgd^ unknown parameters in 
dgd^ equations; there are at most d^L independent so- 
lutions. By solving Eqs. (jSSSp . we can derive a set of in- 
dependent Hermitian operators corresponding to the set 
of their independent solutions. As we have mentioned 
above, for all such Hermitian operators H, its trans- 
pose -ff'^ can be a legitimate Hamiltonian Hse describing 

the dynamics of our triple, (^cIe , | T^Ai ) "8 | T ea2 ) j Hse^ ■ 
Hence, the mission completed. :-) 



IV. A PROTOCOL FOR STATE-STEERING 
TOWARDS THE FULFILLMENT OF THE ME 
CONDITION 

As we have seen in Sec. [TTl as long as our access is 
limited to the principal and ancillary systems S and A, 
the best we can do (from quantum control perspective) 
is to determine an equivalence class of triples in the form 
of (dfi, sea) , Hse)- The information on Hse thereby 
obtained is sufficient to exploit environment as a re- 
source for quantum computing by actively controlling it 
via system S. In order for this Hamiltonian identifica- 
tion to work out, the state \^sea) has to satisfy the ME 
condition in our scenario. 

Therefore, we need a method to steer any given state 
on Hs ^ T~Le towards such a state that fulfills the ME 



condition. In this section, we present a protocol for this 
task allowing us to append an extra (ancillary) system 
Ha- 

As we have mentioned in Sec [H it can be taken for 
granted that the whole state on Hs ®He pure at the 
beginning. In addition, this initial state on Hs ® He 
is assumed to be the same (fixed), but perhaps an un- 
known, state after equilibration. We set t = when each 
iteration of the protocol starts. The protocol proceeds by 
iterating a block of steps that consist of four major ele- 
ments, namely the SWAP operation Q between S and (a 
subsystem of) the time evolution driven by the Hamil- 
tonian Hse, the local filtering operation J^f^p on A, and 
state tomography on S and A. The number of iterations 
of the block is indexed by C. 

Since we need to perform state tomography at suffi- 
ciently high frequency, which we will specify later, and 
there are non-deterministic operations J^f^p, the protocol 
needs to be iterated (ideally infinitely) many times by 
resetting the entire state and the clock. 

A quantity lS.Es a, which plays a central role in de- 
signing the protocol, is defined as a functional of state 
on Hs ® Ha as 



/^EsA := S{psa) - S{pa) + logds, 



(S59) 



where S(p) = — Trplnp is the von Neumann entropy of 
state p [l5l - [l7| . 

The local filtering operation J^f^p is given as J^f^p (p) = 

^lfP^lf, where Flf := -y/Vin~P^, where p']^^ and 
Amin are the inverse of pA on its support and the smallest 
eigenvalue of pA- The local filtering operation succeeds 
with probability Amin • rankp^i. When the local filtering 
fails, we abort the present protocol and restart it from 
the beginning. 

The protocol proceeds as follows: 

Step Q. Ki t — 0, i.e., before any iteration of the fol- 
lowing steps, there is ancillary system A, thus 
dim 'Ha — 0. Set the counter C — Q. 

Step 1. Prepare a standard MES on a pair of new ^5- 
dimensional ancillary systems, Hai ®Ha2- Apply 
SWAPsAi on Hs ® Hai and relabel the resulting 
group of AA1A2 as a new A system. 

Step 2. Apply the local filtering operation J^f^p on Ha 
and increase C by one and call this time tc- 



Step 3. Evaluate AEsa in Eq. (|S59p through state to- 
mography on SA and define ec as 

ec ■-= ^sup {AEsA {t) \ te [tc, tc + At]} (S60) 

during [tc,tc + Atc]- 

Step 4. Terminate the protocol if ec — 0; otherwise, let 
the SE system evolve until AEgA becomes larger 
than (or equal to) ec and go back to Step 1. 
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As mentioned earlier, because state tomography is in- 
volved in Step 3, the evaluation of \Esa can be achieved 
by iterating all the preceding steps many times. Having 
obtained /S.Esa as a function of time over [tc,tc + At], 
we set cc to be a threshold that can be attained within 
this time period. The factor of 1/2 in Eq. (|S60|) is chosen 
merely for convenience to define an achievable threshold. 

How long can At be? It is the time length, within 
which we obtain lS.EgA as a function of time by repeating 
state tomography. As long as our active controls, such 
as SWAP and -^ffp, can be performed in an error- free 
manner, as we assume in this study. At can be arbitrary. 
In order to minimise the overall time length, the shorter 
the A< is the better; however, if At is too short the local 
filtering would succeed with only a very small probability. 
Therefore, a more realistic length of At would be the 
period of 'oscillation' of lS.Es a, which means that the 
value of At may vary from time to time. Nevertheless, 
this strategy should work, since At only needs to have an 
appropriate ec, and its precise, or best, value does not 
have to be known a priori. 

An important observation follows as a theorem: 

Theorem 3. If ds '■= dimHe is finite, the protocol halts 
when C — cIe at latest, i.e., before t = dE^t. In other 
words, there exists a natural number K < ds such that 
AEsEit) — for all t e [tK,tK + At], or equivalently, 
eK = 0- 

Let us first see that AEse at time t is equal to the 
amount of the increment of entanglement between Ha 
and Hs ®He, which is induced when Step 1 of the pro- 
tocol is performed at t. 

Lemma 7. Suppose that the joint system SEA at time 
t in the protocol is described by a state \^sea), o-nd 
I'^'sea) o.'^'d I'^sea) arestates onl-Ls®'HE®'HA®'HAi® 
'Ha2 given as 

\'^'sea) ■■= \^sea)®\Tma.). 
m^A) := SWAPsAA^tEA) 

= \^aa^e)®\Tsa,). (S61) 

where \T A1A2) 0,'^d \TsA2) o-fe written in the standard 
form of ds- dimensional maximally entangled states on 
Hai^Ha2 andy.s<S)HA2, respectively. Then, for psA ■= 
T^eP il'f sea)), 

AEse = £^(1*^^^)) " mn^EA), (S62) 

where EQ"^)) is the amount of entanglement with respect 
to the par tition between Hs ^He and Ha Hai ® 'Ha2 

(Proof) 



mn^iA)) - eh^^'sea)) 

^E{\-^aa^e) «) |TsA.» - E{\-^sea) ® |TaiA.)) 

= «) \TSA2)) - S(|*S£a)) 

= S{pAAi <8) Pmix) - S{pa) 

^^(psa) - S{pa) + logds 
^AEsE- 

In the RHS of the fourth equation, pmix is a completely 
mixed state on Hs- B 
Note that AEse is evaluated with respect to a state 
before the SWAP between S and Ai. During Step 3, the 
state pa ■= T^h^seP il'^ sea)) does not change, and stays 
as a projector on Ha (up to a proportionality constant) 
as a result otJ'{}p. Thus |*1;"e^) is a MES on their sup- 
port, and if AEse > 0, the Schmidt rank of \'S°s1sa) 
strictly greater than that of I'^^s^ea)- This means that 
rankp^T* ~ rankp^ > 1; namely, applying the SWAPsai 
in the Step 1 increases the rank of the state on A by at 
least 1. 

Theorem [3] can now be proven with these lemmas and 
facts. 

(Proof of Theorem [3]) 

Note that the state on the joint system SEA is always 
pure throughout the protocol as long as all the local 
filtering operations succeed. Since the local filtering 
operation J-^p preserves the Schmidt rank of \'^sea) 
with respect to the partition A\SE, the execution of 
Step 1 increases the Schmidt rank of \'^sea) by at 
least 1 (as we have seen above). The Schmidt rank of 
I'^SEa) with respect to the same partition is obviously 
smaller than dE- Thus C is no greater than and the 
statement of the theorem follows. ■ 

The following lemma shows another, more intuitive, 
meaning of AEse- 

Lemma 8. The following equation holds for AEse- 

AEsE = D{pse\\PS ® Pe) + D{ps\\prmx), (S63) 

where D{p\\a) is the relative entropy defined as 
D{p\\a)--^ Trp log p - Trp log a 0, Q, \B, 0, M, M]- 

(Proof) 

AEse 

^S{pe) - S{pse) + logds 

^S{pe) + S{ps) - S{pse) - S{ps) - Trps ^ogPmix 

=IpiS\\E) + D{ps\p,mx) 

=D{pse\\ps Pe) + D{ps\\p7nix), 

where Ip{S\\E)) := S{ps) + S{pe)-S{pse) is the mutual 
information between S and E, and we use the formula 
Ip{S\\E) = D{pse\\ps ® Pe) [3] in the fourth equation. 
■ 

This lemma guarantees that when AEse is small, psE 
is close to ps (81 Pe and ps is almost completely mixed. 
We can also show the following: 
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Lemma 9. When psE = Ps ® pE, there exists a decom- 
position of system A into two distinct subsystems Ai and 
A2; that is, Ha = T~Lai ®'Ha2- Accordingly, a (pure) 
state \^sea) can be written as a product of pure states 
\^SAi) & Hs ® "Hai and I^bAs) G 'H-e ® 'HA2 such that 

I^^sea) = \'^SA,) <S)\^EA2), (S64) 

where \^sAi) and \^eA2) satisfies Tta^P {\^sAi)) = 
ps and Tta2P i\^EA2)) — Pe, respectively. In par- 
ticular, when ps — Prnix, |$SAi) can be chosen as 
a standard maximally entangled state, i.e., \TgAi) — 

dt'^T.tiVs)VA,). 
(Proof) 

Consider any decomposition of A into subspaces A'^ and 
A'2 such that T-La — T~La'-^ <8)'Ha^ and dim'H/i'j^ = dimHs- 
Then there exist the purifications |*i>Ai£;) ^'^'^ \^'eA2) 
Ps and Pe in each subspace, i.e., Tr^i-P (|$5Ai)) ~ Ps 
and TrA^f (I^BAs)) = Pe, respectively. As \'^sea) and 
I'&A^^s) |*&a^e) have the same reduced density matrix 
on T-Ls ® Hsjthere exists a unitary operator Ua on "Ha 
such that 0, 11] 

\^sea) = Ua\^a;s)®Wa'2e)- (S65) 

Defining "Hai ®'Ha2 '■= UaV-a' ®'Ha' proves the lemma. 
■ 

With TheoremOand Lemmas [5] and [51 we finally arrive 
at the following theorem: 

Theorem 4. Let K be the index of the counter C 
when the protocol halts. Then, for all too ^ + ^t, 
{dE,Y^ SEA{tK)) T Hse) satisfies the ME condition for 
[txTtoo), where \^ SEA{tK)) is the state on SEA after 
the K-th (successful) application of the local filtering op- 
eration ^n Step 2 at time tK • 

(Proof ) 

By the construction of the protocol, AEse {t) — for 
all t e [tK,tK + At]. By Lemma [SI psE{t) = ps(t) <E) 



PE{t) and ps{t) is a maximally mixed state for all t € 
[tx, tk + Ai]. Then, due to Lemma[ni for all t & [tK,tK + 
/S.t], there exist Hilbert spaces 7^Ai(0 and HAiit) such 
that Ha ~ 'HAi(t) ®'HA2{t) and \^sEA{t)) satisfies 

TyEP{\^SEAm = P (|TAi(t)5» ® PA2{t), (S66) 

where Ai(t)s) is a maximally entangled state on 
'Hs ®T-LAi{t). Moreover, PA2(t){t) is a projector, because 
of the local filtering operation in Step 2. Therefore, 
{dE,\'^SEA{tK)),HsE) fulfills the ME condition for 
[tK, tK + Ai]. Corollary [2| then lets us finish the proof of 
the theorem. ■ 

We have seen that steering the entire system of A, S, 
and E towards one that satisfies the ME condition can be 
achieved within a finite time (Theorems [3| and [U . Once 
this state-steering has been done, the identification of the 
equivalence class by probing only systems S and A can be 
performed (Sec. IIIip . Therefore, we can always determine 
the equivalence class of a triple associated with a finite- 
dimensional environment, which is essential in identifying 
the Hamiltonian Hse, starting from an arbitrary state on 
the system SEA. 

In the present study, we have assumed that the state 
tomography on SA can be performed perfectly. The fea- 
sibility of the presented protocol for state-steering de- 
pends on this assumption: in order to complete the task, 
we need to check whether AEsE{t) is exactly for all 
t G [tc,tc + At]. Although we have focused on the the- 
oretical aspect of our tomographic scheme for a 'not di- 
rectly touchable system', which is quite remarkable in 
its own right, any quantum operation, including state 
tomography, is always fraught with the effect of unpre- 
dictable noise in reality. Therefore, from a pragmatic 
point of view, we need to modify the protocol and eval- 
uate errors that may occur in the equivalent class iden- 
tification. Since the analysis of errors in the protocol is 
beyond the scope of this paper, we leave it as a future 
project [2^ . 
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